Quantum jumps in a two-level atom 
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A strongly-driven (Q 3> 7) two level atom relaxes towards an equilibrium state p which is almost 
completely mixed. One interpretation of this state is that it represents an ensemble average, and 
that an individual atom is at any time in one of the eigenstates of p. The theory of Teich and Mahler 
[Phys. Rev. A 45, 3300 (1992)] makes this interpretation concrete, with an individual atom jumping 
stochastically between the two eigenstates when a photon is emitted. The dressed atom theory is 
also supposed to describe the quantum jumps of an individual atom due to photo-emissions. But 
the two pictures are contradictory because the dressed states of the atom are almost orthogonal to 
the eigenstates of p. In this paper we investigate three ways of measuring the field radiated by the 
atom, which attempt to reproduce the simple quantum jump dynamics of the dressed state or Teich 
and Mahler models. These are: spectral detection (using optical filters), two-state jumps (using 
adaptive homodyne detection) and orthogonal jumps (another adaptive homodyne scheme). We 
find that the three schemes closely mimic the jumps of the dressed state model, with errors of order 
(7/T2) 2 / 3 , i(7/fi) 2 , and |(7/fl) 2 respectively. The significance of this result to the program of 



environmentally-induced superselection is discussed. 
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I. INTRODUCTION 

The quantum jump, the effectively instantaneous tran- 
sition of an atom from one state to another, was the first 
form of nontrivial quantum dynamics to be postulated 
|EJ. Of course Bohr's theory did not survive the quan- 
tum revolution of the 1920s. In particular, the idea of 
jumps appeared to be in sharp conflict with the continu- 
ity of Schrodinger's wave mechanics In the aftermath 
of the revolution, quantum jumps were revived || with a 
new interpretation as state reduction caused by measure- 
ment. But Wigner and Weisskopf [Q had already derived 
the exponential decay of spontaneous emission from the 
coupling of the atomic dipole to the continuum of electro- 
magnetic field modes. That is, they did not require the 
hypothesis of quantum jumps. Later, more sophisticated 
theoretical techniques, such as the master equation, were 
developed for dealing with the irreversible dynamics of 
such open quantum systems [f>|-|7|. In the master equa- 
tion description, the atom's state evolves smoothly and 
deterministically. Perhaps as a consequence, interest in 
quantum jumps as a way of describing of atomic dynam- 
ics faded. 

Quantum jump models for atoms were never entirely 
forgotten; the dressed state model || was used success- 
fully to give an intuitive explanation of the Mollow triplet 
H in resonance fluorescence. However, it was the elec- 
tron shelving experiments of Itano and co-workers |lC| ] 
which refocussed attention on the conditional dynamics 
of individual atoms. Subsequent work on waiting time 
distributions 0,021 led to a renewal of interest in quan- 



tum jump descriptions |13|. It was shown by Carmichael 
|l4| that quantum jumps are an implicit part of stan- 
dard photodetection theory. This link between continu- 
ous quantum measurement theory and stochastic quan- 
tum evolution for the pure state of the system was consid- 
ered by many other workers around the same time and 
subsequently |l5| "|2^]. Independently, Dalibard, Castin 
and M0lmer [|25|| derived the same stochastic Schrodinger 
equations, driven by the need for efficient methods for nu- 
merically simulating moderately large quantum systems. 
This technique, called Monte-Carlo wavefunction simula- 
tions, has been applied to great advantage in describing 
the optical-cooling of a fluorescent atom [^6-30|. Re- 
gardless of the motivation for their use, the evolution of 
systems undergoing quantum jumps (and other stochas- 
tic quantum processes), are known widely as quantum 
trajectories Q| . 

Remarkably, also around the same time as quantum 
trajectory theory was being developed, an entirely dif- 
ferent theory of quantum jumps in atomic systems was 
proposed by Teich and Mahler (TM) These au- 

thors claimed that the quantum jumps in their theory 
also corresponded to photon detections. However, it is 
plain that, except in trivial cases, the TM trajectories 
are different from the quantum trajectories from direct 
photodetection. The internal consistency of the TM tra- 
jectories has also been criticised by one of us but 
it turns out that there are some subtle issues involved 
here (as will be discussed later). This, combined with 
the simplicity of the TM theory, suggests that it may be 
worth a closer look. 



*Electronic address: wiseman@physics.uq.edu.au 



1 



In this paper we use the TM trajectories as a base 
from which to explore quantum jumps for the simplest 
non-trivial atomic system, a two-level atom in a strong 
resonant driving field. First (Sec. II) we summarize the 
TM theory and, following Teich and Mahler, derive the 
TM trajectories for this example. We then identify vari- 
ous features of these trajectories, and the interpretation 
given them by Teich and Mahler, and ask the question: 
can any of these features be reproduced by other jump 
models? Specifically, we examine four other jump mod- 
els: the dressed atom model (Sec. Ill), spectral detec- 
tion (Sec. IV), two-state jumps (Sec. V), and orthogo- 
nal jumps (Sec. VI). The last three models derive from 
quantum trajectories but use detection schemes which 
are progressively harder to implement. The results of 
these comparisons are discussed in Sec. VII. 



II. THE THEORY OF TEICH AND MAHLER 

A. General Theory 

In this section we summarize the theory of Teich and 
Mahler in our own notation. The most general form of 
Markovian master equation for a quantum system is the 
Lindblad form [02l 



f)= Cp = £ 1C vP + ArrP- 

Here the reversible and irreversible terms are 



CicvP = -i[H, p], 



(2.1) 

(2.2) 
(2.3) 



Here H is an Hermitian operator, the Cj are arbitrary 
operators and I? is a superoperator defined for arbitrary 
operators A and B as 



V[A]B = ABA^ - \{A^A, B}. 



(2.4) 



Say the solution of the master equation Eq. ( |2.l| ) is 
p(t). Then p(t) can be diagonalized at any time as 



p(*) = $>p(*)I0i.(*)Wp(*)I, 



(2.5) 



where D is the dimension of the Hilbert space of the sys- 
tem. The time-dependent eigenstates of p(t) are always 
orthogonal to one another, and so evolve according to a 
Unitary transformation: 



0_ 

Of 



\4>Jt)) = -iK(t)\<f>Jt)), 



(2.6) 



for some Hermitian operator K(t) which will depend 
upon p(t) but not upon //. The rate of change of the 
eigenvalues is thus given by 



(2.7) 



= -iTr{p(t)[K(t),\ ( f >ll (t))(^(t)\}} 
-zTr{\^(t))(^(t)\{H,p(t)}} 

+Tr[\^(t))(^(t)\C irT p]. (2.8) 

Using the cyclic properties of the trace operation and 
the fact that [|(/) At (i))(0 M (i)|, p(t)] = 0, the first two terms 
vanish. Thus one is left with 

d 



ic 1 

2 U J 



c)c ]P {t)-\p{t)c]c 3 |<^(i)>, (2.9) 



where the Einstein summation convention for latin in- 
dices is being used. Using Eq. ( |2.5| ) and the completeness 
relation 1 = J2„ \^v{t)){^>v{t)\, one obtains finally 

J^vw = E [<^(*)|c5i^(*)><^(*)|c,-i^(*)>^(*) 

V 

- (^t)\c]\Mt))(Mt)\cj\Mt))p^)} ■ (2-io) 



Teich and Mahler interpret the state matrix p(t) as 
pertaining to an ensemble of individual systems. The in- 
dividual systems, they say, are always in one of the states 
\4>fi{t)) but jump stochastically between these state s. To 
derive the rates of these jumps, they rewrite Eq. ( 2.1C ) 
as 



where 



iW*)=Ei<^(*)i*(*)>r 



(2.11) 



(2.12) 



In this form it is clear that R^ v can be interpreted as the 
probability per unit time for the system in state |0„(i)) 
to jump to state |</> M (£)). Note that jumps which leave 
the system in the same state also occur, since is in 
general non zero. Between jumps, the system remains 
in one of the states \cj> u (t)) , an d hence changes smoothly 
in time according to Eq. ( |2.6| ). It is clear that the TM 
traject ories do reproduce the mast er e quation Eq. fl2.ll) , 
as Eq. (2.11) is equivalen t to E q. (2.1). 

From the form of Eq. ( |2.12 ) it appears that the rates 
may depend on the way that C lrv is written. The choice 
of operators Cj in the definition of C lrT is not unique; a 
unitary rearrangement leaves invariant. That is, if 
we define 



c 'k — UkjCj 



where UkjUt, then 



(2.13) 



(2.14) 
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It turns out that this rearrangement leaves un- 
changed also: 



k 



1(^141^)1 > 



(2.15) 



so that the TM trajectories do not depend on the way 
that £j rr is written. 

However, it turns out that the rates R^ do depend on 
the way that C is split into C rcv and C m . It is possible to 
change £ rev and £; rr while leaving L unchanged by the 
following transformation: 



cj + Xj, 



H -> H' 



H - -[\* Cj - \ jC ]} 



(2.16) 
(2.17) 



where the Xj are arbitrary c-numbers. It is simple to 
see R^ v is not invariant under this transformation. Tc- 
ich and Mahler do not state what A., should be chosen 
before applying their method; their Eq. (1) already as- 
sumes the separation C into £ rov and C llr . One general 
choice which one might make (and which coincides with 
the choices made by Teich and Mahler in the simple sys- 
tems they consider) is that the Xj be such that c'j be 
traceless operators. 

In steady state, p(t) is time-independent and we will 
denote it simply 



(2.18) 



In this case the rates are also time-independent and are 
given by 



RlLV 
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(2.19) 



In what follows we will only consider this stationary 
stochastic evolution. 



P = 



ft 2 + Qja v + 7 2 (1- a z )/2 
2ft 2 + 7 2 



(2.21) 



where the Pauli matrices have their usual meaning. To 
simplify matters, consider the strong driving limit ft 3> 
7. Then to first order in 7/ft we have 



ft + jo-y 

2ft 



(2.22) 



Diagonalizing this p as in Eq. (2.18|) yields the following 



V2|0i) = |s>-<|e>, 
y/2\<h) = \g)+i\e), 



1 



P2 



1 



Following the choice in Eq. (2.20) of £; r 
jump rates can be written 

R fiU = ~/\{<t>M<t>»)\ 2 - 



(2.23) 
(2.24) 

(2.25) 
(2.26) 

jV[a], the 
(2.27) 



Substituting the expressions for \4>i^) into Eq. (2.27) 
gives 



R11 = -R12 — R21 = Ri 



7/4. 



(2.28) 



The fact that R21 = R12 means that the TM trajectories 
would actually predict pi = f>2 = 1/2. To obtain rates 
which would give the correct probabilities it would be 
necessary to determine |0i 2) to higher order in 7/ft. 

Teich and Mahler associate each jump, including those 
which leave the system state unchanged, with the emis- 
sion of a photon. In this case the total rate of photoe- 
mission is therefore 7/2. To first order in 7/Q this agrees 
with the expression from the stationary density operator: 



B. The Two-Level Atom 

Consider an atom with two relevant levels {\g}, |e)}. 
Let there be a dipole moment between these levels so 
that the coupling to the continuum of electromagnetic 
field modes in the vacuum state will cause the atom to de- 
cay at rate 7. So that the atom does not simply decay to 
the state \g), add driving by a classical field (such as that 
produced by a laser) of Rabi frequency ft. Then the evo- 
lution of the atom's state matrix can then be described 
by the master equation (in the interaction picture) 



,ft r 



, p] + ~fV[a]p. 



(2.20) 



Here a = \g)(e\ and o~ x = a + . 

The stationary stationary solution of Eq. (2.20) is 



7(e|p|e) 



2 
2 



o 



(2.29) 



Photons associated with state-preserving jumps Teich 
and Mahler assign to the central peak of the Mollow 
triplet, while the state-changing emissions they assign to 
the sidebands. The ratio of intensities in the central peak 
and the sidebands thus agree with those in the Mollow 
triplet §. 

The TM trajectories have a number of interesting char- 
acteristics. For the purposes of the rest of the paper, 
there are three in particular to which we wish to draw 
attention: (i) In steady state, the jumps supposedly cor- 
respond to photoemissions into the three peaks of the 
Mollow triplet; (ii) In steady state, the atomic state is 
always one of two fixed pure states; (iii) At all times, the 
state after a jump is orthogonal to the one before. 
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III. THE DRESSED ATOM MODEL 

The dressed atom model [0 can be applied to an ar- 
bitrary atomic system but for simplicity we will con- 
sider only the case at hand, the resonant two-level atom. 
The mod el is based on replacing the Hamiltonian in 
Eq. ( |2.2C| ), in which the driving laser is treated as a clas- 
sical field, with one in which the driving laser is treated 
as a single-mode quantum field. Putting in the self- 
Hamiltonians of the atom and field gives (in the rotating- 
wave approximation) 



^(aV + aat), (3.1) 



where a is the annihilation operator for the driving field 
and g is the dipole coupling constant also known as the 
one-photon Rabi frequency. This Hamiltonian has eigen- 
states 



v / 2|n,±) = |n)| 5 )±|n-l)|e), 



(3.2) 



where \n) are number states of the driving field. These 
eigenstates are known as dressed states of the atom, as 
opposed to {|e), \g}}, the bare atomic energy states. They 
have energies 



E n< ± = h(nujQ ± \/ng/2). 



(3.3) 



For large coherent driving field and small coupling con- 
stant g the classical approximation is valid and one can 
replace \/ng by \ffig = tt . Then, for n ~ n, the ladder 
of energy eigenstates fl3.3j ) will consist of pairs of closely- 
spaced rungs, with an inter-pair separation of Hujq, and 
an intra-pair separation of Ml. 

Now one can interpret the Mollow triplet in terms of 
spontaneous-emission-induced transitions between these 
stationary states. If the dressed atom is in one of the 
states |n, ±), it can spontaneously emit a photon and 
drop down a rung on the ladder. If it drops to \n — 1, ±) 
(that is, the atom effectively remains in the same state), 
then the change in energy of the dressed atom is fkoo and 
so the frequency of the emitted photon must be u>q — in 
the central peak of the triplet. If the atom changes state 
via a transition to \n — 1,^)) then the frequency of the 
emitted photon must be loq ± £1 — in the sideband peaks. 

The rates of these transitions are calculated according 
to the same formula as in Teich and Mahler's scheme. 
For example, if the system were in the dressed state |n, +) 
then the probability of it to spontaneously emit and make 
a transition to the dressed state |n — 1, — ) is 



Substituting this expression into Eq. (3.2) 
similarly for the other transitions, gives 



R. 



= R++ = = R++ = 7 /4. 



(3.4) 
and doing 



(3.5) 



That is, all of the transition rates are equal. A strict or- 
dering of transitions will also occur, ensuring that a high 



frequency photon must be emitted between emissions of 
low frequency photons, and vice versa. 

All of the features of the dressed-state model discussed 
above agree with those of the TM theory. Teich and 
Mahler even go so far as to call their theory a "general- 
ization of the dressed state picture." |3l[] However, in one 
crucial respect the two theories are completely at odds: 
the state of the atom. In TM t heory it ju mps between 
the states and \<j> 2 ) of Eqs. fl2.23p and ( |2.24| ). These 
are eigenstates of a v . By contrast, for a coherent driving 
field the dressed states of the atoms are 



-M 2 /2. 



(|n)| 5 )±|n-l)|e)) 



(3.6) 



For \a\ large the field and atom states are very nearly 
not entangled, and the state of the atom is very close to 
a pure state |±). For a real, as required to make the 
replacement 



a, + 

2 (aV 



acr T ) — (ct + ct 1 ), 



(3.7) 



so as to yield the master equation ( 2.20| ) with £1 = ga, 
these pure states are defined by 



V2\±) = \g}±\e) 



(3.8) 



These are eigenstates of a x . Thus the atomic states in 
the dressed atom model are as different as it possible to 
be (on the Bloch sphere) from the atomic states in the 
TM model. 



IV. SPECTRAL DETECTION 



The disagreement between the dressed state model and 
the TM model, both of which claim to describe emis- 
sion into the three peaks of the Mollow triplet, suggest 
that it would be worthwhile investigating such emissions 
by a third, independent theory. The theory we use in 
this section is that of quantum trajectories, developed 
initially by Carmichael jl4| and subsequently by many 
other authors. This theory is essentially an application of 
quantum measurement theory of continuously monitored 
systems, so we will briefly review this theory. We will 
consider only efficient measurements in which no infor- 
mation is lost. In the optical context, this requires com- 
plete collection of the emitted light and unit efficiency 
photodetectors. 
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A. Quantum Trajectories 

The aim of quantum measurement theory is, given the 
initial state of the system, to be able to specify the prob- 
ability of a particular measurement result, and the state 
of the system immediately following this result. Say the 
measurement result is a, a random variable which will be 
assumed discrete for convenience. Then both the proba- 
bility and the conditioned state can be found from the set 
of operators {U a (T)}. Here T is the duration of the mea- 
surement and the operators tt a (T) are arbitrary, with one 
condition, 



£nt(T)n a (T) = i, 



(4.1) 



where the sum is over all possible results. This is known 
as the completeness condition [ |33|| . 

The probability for obtaining a particular result a is 
found from the measurement operator by 



where 



P a = Tr[p a {t + T)], 



p a (t + T) = n a (T)p{t) fit(T) 



(4.2) 



(4.3) 



is an unnormalized density operator, where p(t) is the 
density operator at the beginning of the measurement. 
The state of the system conditioned on the result a is 
simply given by 



p a (t + T) = p a (t + T)/P a . 



(4.4) 



This economy of theory is a consequence of a more fun- 
damental notion of probability relating to projectors in 



Hilbert space |34 
pure (p = \ip)($ 
state is obviously 



If the initial state of the system is 
then the unnormalized conditioned 



\4> a (t + T)) = n a (T)\ip(t)). 



(4.5) 



If the measurement is performed, but the result ignored, 
then the new state of the system will be mixed in general 
and cannot be represented by a state vector. It is equa l 
to the sum of the conditioned density operators (4.4) 
weighted by the probabilities (|4 



(4.6) 



p(t + T) = ^P Q p Q (t + T) 

a 

= Y / n a (T)p(t)ni(T). 



It is easy to verify from the completeness condition (4.1) 
that Tr[p(i + T)] = 1, as required by conservation of 
probability. 

Continuous measurement theory can now be cast as 
a special case of quantum measurement theory where a 
constant measurement interaction allows successive mea- 
surements, the duration of each being infinitesimal. If 



the state matrix at time t is p{t), then the unnormalized 
conditioned density operator after the measurement in 
the interval [t, t + dt) is denoted 



t (t + dt) = n a (dt)p(t) ni(dt), 



(4.7) 



where the operators Q a (dt) are arbitrary as yet. The un- 
conditioned infinitesimally evolved state matrix is then 



p(t + dt) = J2 tt*{dt)p{t) n+ (dt). 



(4.8) 



This represents the evolution of the system, ignor- 
ing the measurement results. If the fl a (dt) are time- 
independent, this nonselective evolution is obviously 
Markovian (depending only on the state of the system 
at the start of the inter val) . 

As noted in Sec. [I A , the most general form of Marko- 
vian equation of motion for the density operator of a 
system is a master equation o f th e Lindblad form [j32| . 
Consider a special case of Eq. (2.1) where there is a sin- 
gle Lindblad operator c. Then, by inspection, there are 
only two necessary measurement results (say and 1), 
with corresponding operators 



(4.9) 
(4.10) 



D,i(dt) = Vdtc, 

n (dt) = 1 - [iH+^c] dt. 



The two unnormalized conditioned density operators are 

px(t + dt) =dtcpc f , (4.11) 
p {t + dt) = p + dt (-i[H, p] - i{c + c, p}) . (4.12) 

It is easy to see that 

p(t + dt) = fa{t + dt) + p {t + dt) = p + dtp, (4.13) 



where p is given by the master equation (2.1). 

Evidently, almost all infinitesimal intervals yield the 
measurement result 0. Upon such result, the system state 
evolves infinitesimally (but not unitarily) . Whenever the 
result 1 is obtained, however, the system state changes 
by a finite operation. This discontinuous change can be 
justifiably called a quantum jump, and the measurement 
event a detection. As with all efficient measurements, 
if the initial state of the system is pure, then the condi- 
tioned state of the system will remain pure. The stochas- 
tic evolution of such a conditioned state has been called 
a quantum trajectory by Carmichael [|l4| , as discussed in 
the introduction. 

Now consider the application of continuous quantum 
measurement theory to the damped, driven two-level 
atom which is the subject of this paper. The master 



equation for this system is given by Eq. (2.2C), and this 
can be most simply unraveled using the two measurement 
operators 



f2o(<it) = 1 — i—o~ x dt — — a'adt, 
fli(dt) — \J jdt a. 



(4.14) 
(4.15) 
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In this unraveling, the rate of jumps is 7(0^(7), which is 
identical to the rate of spontaneous emissions. Thus this 
unraveling corresponds to direct detection (by a photode- 
tector) of all of the light emitted by the atom. Indeed, 
the detection operator £l\(dt) is proportional to the part 
of the field radiated by the atom. 

It is evident that under this detection scheme, the 
atomic state immediately following a jump is always the 
ground state \g). Thus these quantum jumps are quite 
different from those in either the dressed state or the 
TM model. This is not surprising, as the photodetector 
does not distinguish photons of different frequencies, as 
required in these other models. 

In order to distinguish photons of different frequency 
it is necessary to use optical filters in the detection ap- 
paratus prior to photodetection. There are a number of 
ways to describe such a detection scheme. One is to use 
a non-Markovian evolution equation for the system, as 
done by Jack, Collett and Walls p5| . However, this has 
the added complication that the system is represented 
not by a single state vector but by a sequence of state 
vectors (which in principle is infinite). Also, this for- 
malism does not directly give the state of the system at 
the time t of the measurement, but rather the retrodicted 
state 36 1 at time t — T m . Here T m is the time required 
by the filtering process. In this work we adopt a different 
procedure in which the filters are described as quantum 
optical systems (cavities). In our formalism only a sin- 
gle state vector is required, but it exists in the enlarged 
Hilbert space of the system (the two-level atom) plus fil- 
ters. This state vector applies to the system plus filters 
at the moment of measurement (ignoring the propagation 
time for light between the various elements). Finally, our 
method is amenable to an approximate analytic solution 
in a suitable limit. 



B. Cascaded quantum systems 

In order to describe auxiliary quantum systems (fil- 
ters) as part of the detection scheme it is necessary to 
use "cascaded systems theory" as it has been called by 
Carmichael Cascaded systems are different from 

coupled systems, because the interaction only goes one 
way. That is to say, the first system influences the sec- 
ond, but not vice versa. One mechanism for achieving the 
required unidircctionality is the Faraday isolator which 
utilizes Faraday rotation and polarization-sensitive beam 
splitters. This is most practical for the case of cavities, 
where the output is a beam of light. In principle this 
technique could be applied to the radiation of an atom as 
well, if the atom were placed at the focus of a parabolic 
mirror so as to produce an output beam, as shown in 
Fig. I. 

A quantum theoretical treatment which incorporates 
this spatial symmetry breaking at the level of the Hamil- 
tonian has been given by Gardiner [R7| . If the propaga- 



tion time between the source system the and driven sys- 
tem is negligible then a master equation for both systems 
may be derived. This result was obtained simultaneously 
by Carmichael (|T), who used quantum trajectories to il- 
lustrate the nature of the process. Since we wish to de- 
scribe the monitoring of the outputs of the filter, we will 
follow Carmichael's approach. 

Let the first system be a two level atom obeying the 
usual master equation 



-i[(n/2)<T x ,p]+V[yfy<r]p. 



(4.16) 



As noted above, the field radiated by this atom is rep- 
resented by the operator y/jcr. Now say this field (plus 
the accompanying electromagnetic vacuum fluctuations) 
impinges upon the front mirror of a Fabry-Perot etalon 
of linewidth 2T. That is to say, the intensity decay rates 
through the front and rear mirror are both T. Then the 
master equation for the state matrix W for the combined 
system is 



W 



(tt/2)a x + UaC^a + i^/jf(a f a - era t ), W 
V[^o + Vf a]W + V[Vr a]W, (4.17) 



where u a is the detuning of the relevant mode of the 
etalon (relative to the atom and its resonant driving field) 
and a is its annihilation oper ator. 

It can be verified from Eq. ([4.17 ) that tracing over the 



cavity mode a yields Eq. ( l.lq ) for the atom alone. That 
is to say, the filter does not directly affect the atom, which 
is as desired. The apparent coupling term in the Hamil- 
tonian in Eq. (4.17) is canceled by the interference in the 
irreversible term with Lindblad operator y/j a + \/T a. 
This operator represents the radiated field from the front 
of the resonator. It is the Faraday isolator or equivalent 
mechanism which prevents the interaction of this field 
with the atom. The second Lindblad operator VTa rep- 
resents the field radiated from the rear of the resonator. 

Since the filter produces two output fields (that passed 
and that rejected), monitoring the system requires two 
photodetectors. Hence the system will now be described 
by three measurement operators, £l\(dt) corresponding 
to the detection of a "rejected" photon (off the front of 
the etalon), tt a (dt) corresponding to a "passed" photon 
(from the rear), and Qo(dt) corresponding to no detec- 
tion of a photon in the interval [t, t+dt). These operators 
are given by 



^(dt) 

n a (dt) 




n 



ilo(dt) = 1 — idt yjj a a) a + ~^ a x 
It is easy to verify that 



dt + Ta^a + y^yTa'a 



(4.18) 
(4.19) 



(4.20) 



G 



J2 £i a (dtfwn a {dt) = wdt, 



(4.21) 



a=0,l,o 



where W is given by Eq. (4.17) 



C. Two-Level Atom with One Filter 

We wish to consider now the specific case where the 
filter cavity is designed so as to pass photons from the 
high frequency peak of the Mollow triplet, while rejecting 
photons from the middle and low frequency peaks. In the 
limit f2 3> 7 (which is required for the three peaks to be 
well separated), the upper peak is centered at frequency 
u>o + il. Thus we choose the detuning of the etalon to 
be uj a — Q. The two measurement operators £l±(dt) and 
S} a (dt) are unchanged, while flo(dt) then becomes 

tt (dt) = 1 - idt^ (a x + 2a) a) 



dt(J-a ] o + Ta*a + VtTVct) . (4.22) 



Now in order to pass almost all of the high frequency 
photons but reject almost all of the middle and low fre- 
quency photons we require a filter bandwidth 2r sat- 
isfying Q, ^> r ^> 7. In this limit the cavity relaxes 
much faster than the rate of spontaneous emission by 
the atom. Since on the cavity time scale, photodetections 
are infrequent events, it makes sense to consider the ba- 
sis which diagonalizes the no-jump operator flo(dt). Be- 
cause this operator is not normal (that is, it does not 
satisfy [fio(cfa)^, fta(dt)] — 0), its eigenstates are not or- 
thogonal. Nevertheless they form a complete basis and 
to zeroth order in (7/fi) and (7/T) they are orthonormal. 
The exact expressions for the eigenstates are 



\s±) = J2( h tj\ h )\j)+ l ij\m)- 

j=n 

Here \h) and \l) are atomic states defined by 
\h) = Ml.9> + v\e), 

\i) = v\g) -He), 



where 



4fi 2 



(4.23) 



(4.24) 
(4.25) 



(4.26) 



Note that these atomic states are not exactly orthogonal. 
However, with an error of ord er (7/T2) 2 they are equal to 
the dressed states of Eq. ( |3.8|) : 



\h) = 

10 = 



o( 7 2 /fn 



(4.27) 
(4.28) 



and hence are very nearly orthogonal. In Eq. (4.23) the 
states |j) are eigenstates states of a^a. The coefficients 
h^j , l^j are defined by the recurrence relations 



nJ h + hj - CTn 



e, = 



T + (p-l)£ J -_i], (4.29) 



nJ h + hj - On 

for j > n where 



P 



and the initial conditions for each recurrence chain are 



fW 



(4.30) 



(4.31) 



h + =1 l + 

,tJ nn ' nn 

form ~ , l nn 



0, 

1. 



In Eqs. (|4.29| ) and ( |f.3C| ) 



0-+ = h + n\ a , 
a~ = Xi + n\ a , 



and 



h 

h = 



4 2V 4 ' 

— r — iuj a = — r — in, 



(4.32) 
(4.33) 



(4.34) 
(4.35) 



(4.36) 
(4.37) 

(4.38) 



As well as appearing in the recurrence relations, the 
define the eigenvalues: 



n (dt)\st) = (i + *tdt)\st). 



(4.39) 



Before proceeding further we make the assumption 
that terms of order (j/fl) 2 may be ignored. This will 
be justified later. Under this assumption /1 = v and 



p = q = i so that 



l^t) = E(^l+)^ + ^l-)b-» 



(4.40) 



where 



"j 



h + hj — On 

h + h j - o",f 



«i-i-C--i)/ 2 ' 



h + hj - On 

where 

7 «fi 

A a = -r - iCl, 



(4.41) 
(4.42) 



(4.43) 
(4.44) 

(4.45) 
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The next assumption we make is that we need consider 
only the two states \Sq). This is based on the observa- 
tion that the real part of is j/A+nT. This means that 
if the state is prepared in a superposition of states |S^), 
the decay rate for the component with n > is much 
faster than for that with n = 0, since r 3> 7. Thus, 
given that no detection occurs, the system will soon find 
itself in the subspace spanned by the states \Sq). We 
will return later to the question of how much error is in- 
troduced by this approximation. Meanwhile, solving the 
recurrence relations ( |4.41 ), (4.42) yields 



V7 



|s + ) = l+)|o) - + (4 - 4(i » 



-)|o)-^Pl + )H)-^l-)ii) 



2Q 



(4.47) 



where the terms ignored with 2 or more photons in the 
cavity of order 7/O. Note that when the system is in 
state \S^) the atom is substantially in the dressed state 
l±>- 

Consider the case when the total system is in state 
\Sq). It will remain in this state until a detection oc- 
curs. If a photon passed through the cavity is detected, 
the new state will be proportional to 



»T, 



n a (dt)\s+)^\-)\o)--\+)\o) + o 



I5'n 



o 

-IS 4 



n 



o 



/7T 



|1) (4.48) 



(4.49) 



That is, the detection of a high-frequency photon trans- 
fers the atomic state from being predominantly in the 
high energy dressed state |+) to being predominantly in 
the low energy dressed state |— ). This is in line with 
the simple dressed state model. Moreover, if instead a 
rejected photon (from the middle or lower peak) is de- 
tected, the new state is proportional to 



tti(dt)\S£) oc |+)|0) +0 



ID 



= ls n H 



(4.50) 
(4.51) 



That is, the atomic state is substantially unchanged. 
This is again as expected from the dressed state theory, 
as a low frequency photon from an atom in the dressed 
state |+) would be impossible, and a resonant frequency 
photon would leave the atomic state unchanged. 

The picture so far from this detection scheme is re- 
markably close to the dressed state model. However, 
the correspondence breaks down when we consider the 
system initially in the state \Sq). Then when a passed 
photon is detected the new state is 



n o (^)|5 -)oc^|+)|0) + |-)|0) + 



|1>, (4.52) 



which is not close to either dressed state. Similarly if a 
rejected photon is detected the new state is 



ni(<ft)|So>oc|5>|o> + o(r/n). 



(4.53) 



which is an equal superposition of the two dressed states. 

The reason for this failure of the dressed state model is 
that the detection scheme does not distinguish between 
photons in the central and lower peak of the Mollow 
triplet. It might be thought that tuning the cavity to the 
central frequency loq could solve this problem, as then a 
passed photon would leave the atom in the same dressed 
state, while a rejected photon would swap the atom from 
one dressed state to the other. This does work for short 
times, as shown in Fig. 2, if the atom starts in a dressed 
state. However, because there is no way of distinguish- 
ing between high and low frequency photons, errors ac- 
cumulate and soon the experimenter could not tell which 
dressed state the atom is in. In theory, the atom is still 
approximately in a pure state, but to know what pure 
state it is in would require timing of the detections to 
time scales less than . This would be difficult in 
practice and is also counter to the spirit of the dressed 
state model. 

To attempt to properly replicate the dynamics of the 
dressed state model (or, perhaps, the TM model) it is 
necessary to distinguish all three peaks of the Mollow 
triplet. This requires two filters and is investigated in 
the following section. 



D. Two-Level Atom with Two Filters 

Consider the set up shown in Fig. 3 with two cascaded 
filter cavities with annihilation operators a and b. The 
master equation for this system is 

W = -i [(0/2)0-3 + cu a a f a + LUbtfb, W] 

— iVT [y/j (ia^a — iaa)), W] 

- iVf [i(y/j(T + Vfa) f b - i(V7cr + Vf a)b f , W] 
+ V[^a + Vfa]W + V[Vfa]W 

+ V[^a + VTa + VTb]W. (4.54) 

Here we have chosen the bandwidth of the second cavity 
to equal that of the first, 2r. 

In this case we have three detection events. If we 
choose uj a — +0 (as before) and ut, = — 0, then photons 
passed by filters a and b will be from the high and low 
sidebands respectively and those rejected by both cavi- 
ties will fall in the central peak. The four measurement 
operators are, in the usual interaction picture, 



Oi(rft) = \J jdt <r 
Q a (dt) = Vfdia, 
Q b (dt) = Vfdib, 



(4.55) 
(4.56) 
(4.57) 
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Q (di) = 1 - idt— (a x + 2a) a - 2b%) 
— dt ( —a* a 



dt i 



Cy/jT a t a + Tb t a + y/jTb^aj . (4.58) 



Note that detection of a rejected photon now involves 
field amplitudes emitted from the atom and both cavi- 
ties. 

To attack the problem we again find the eigenstates of 
Q,o(dt). For n, m natural numbers these are given by 

\sD = EE (h± mjk \h)\m + iLjMm 

j=0 k=0 

(4.59) 

The recurrence relationship for ^ nm „- fc and l^ m j k is 
slightly more complicated than for the single cavity case. 



nmjk 



(X h + Xaj + X b k - (T± m ) 




i-l,k + (P " ^ 



± 

nm,j — l,k 



nmjk 



+ Ty/k(j + l)h± mJ+ltk _ 1 } 
(A; + X a j + X b k - <7±„ 



(4.60) 




(4.61) 



Here A Q = — T — ift as before while At, = — T + til, and 
the eigenvalues of flo(di) are 1 + a^; m (dt) where 



i = A;, + A Q n + X b m, 
CT nm = A/ + A a n + X b m. 



(4.62) 
(4.63) 



The boundary conditions for the recurrence relations are 



h + =11+ =0 

nmnm ' nrnnrn ' 



nmnm 



? ^nmnm 



(4.64) 
(4.65) 



Note that there is an asymmetry in j and fc in Eqs.(4.6C) 
and (4.61) due to the ordering of the cavities. This is also 
manifest in the the range of j and k for a given n and m. 
Starting at k — m, j should range from j = n + 1 to oo, 
while k is then continually incremented, and for each k 
value, j should run from max(n — k + m, 0) to infinity. 

As in the single cavity case, we are interested in the 
longest-lived states \Sq ). Also, we make the the same 
approximations stemming from the limits ft T 7. 
Under these approximations we find the following expres- 
sions (where the 00 subscript has been omitted) 



IS" 1 



-)|00) 



2VF 



-)|10) 



r7, 



2ft 



(|10) - |01)) 



^l+>|ii> + o(g 



o 



7 3/2 

r 3 /2 



(4.66) 



|5-) = |-)|00) + ^L|+)|01) 

+ ^P(2|-)|01)-2|-)|10)-|+)|10)) 



4ft 



sr 1 



l-)UD 



o 



n 



/~3/2 





(4.67) 



Here the omitted states of order 7/ft have two photons 
in one cavity and none in the other, while those of or- 
der (7/r) 3 / 2 have two in one and one in the other. Note 
the asymmetry in the terms of order ^/T^j /ft due to the 
ordering of the cavities. 

Now imagine the total system is in state \S + ). It will 
remain in that state until a detection occurs. The new 
state conditioned on the detection of a photon passed by 
cavity a is 



n a (dt)\s+) oc|-)|oo) + -^|+)|oo) 



4VT 



->|01) + O 



ft 



(4.68) 



Note that to zeroth order, the new system is in state 
|— )|00), as expected: 



Sl a (dt)\S + ) cx \S~ 



(4.69) 



Moreover, the rate for this detection to occur is 

(s + \ni(dt)n a (dt)\s+)/dt, (4.70) 

which to zeroth order evaluates to 7/4, as expected from 
the dressed atom model. 



Returning to the more complete expression (4.68) for 
the system state following an a detection, there is an er- 
ror due to the second term, which is the amplitude for 
the system to jump to the wrong dressed state. The 
magnitude of this error is clearly 



ft 2 ' 



(4.71) 



As well as this error, the correlation between the b cav- 
ity and the atomic state is already present in the third 
term |+)|01), with the same sign as in the entangled state 
IS* - ). Thus any subsequent detection through the b cav- 
ity will have the same effect as if the system were in state 
IS* - ), namely to put the system back in the state |+) |00). 
However, if a photon rejected by both cavities is detected 
before the third term has decayed to its stationary value, 
the system will not jump into the correct appro ximat e 
dressed state. This can be easily verified from Eq. (4.68). 
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Since the rate of such detections is of order 7/4 and the 
square of the amplitude for the third term decays at rate 
2r, the error introduced by this transient scales as 



^transient 



_L 
sr' 



(4.72) 



A different sort of error occurs when a photodetection 
occurs which would be forbidden by the dressed state 
model. In the present situation, when the system starts 
in the state \S + ), the forbidden detection is through cav- 
ity b. Since the (unnormalized) state conditioned on this 
detection is 



n b (dt)\s+) = vW2| 



-iT 
V2Q 



|00) 



V7 



4V2r 



1 10) 



(4.73) 



the probability for this detection to occur scales as 



^forbidden 



r 2 

2ft 2 



7 

32r' 



(4.74) 



Turning now to the other allowed detection, we find 
that to zeroth order 



n 1 (dt)\S + ) oc \S + ), 



(4.75) 



and the rate of this process is again 7/4. Similarly, if the 
atom starts in the state \S~) the to zeroth order 



n b (dt)\s-) oc |5+), 
n^dt^s-} a \s~), 



(4.76) 
(4.77) 



and the rates are 7/4, while the probability of a forbid- 
den d etection through cavity a is similar to the expression 
( 4.74 ). Thus for 7 <C T <C ft, the system almost always 
jumps between the states \S^). This is confirmed by the 
numerical simulations shown in Fig. 4. 

The final source of error is that in these states the 
atomic state is n ot th e pure dr essed state expected. 
Rather, from Eq. ( [4.66 ) and Eq. (4.67), the orthogonal 
dressed state is entangled with the cavity states, with am- 
plitude ^Fj /2VT . Thus the probability for not finding 
the atom in the expected dressed state scales as 



^entangled 



4r 



(4.78) 



There is another error introduced by the fact that 
the atomic states \h), \l) differ from the dressed states 
|+), |— ) by an amount of order (7/ft) 2 . However since 
r>7 this is negligible compared to other errors noted 
above, such as e wrong . 

To determine the total probability for deviation from 
the predictions of the simple dressed state model, we add 
the four sources of error discussed above. The result is 



T 2 7 
Etotal = a--r + p—, 



"ft 2 



4r' 



(4.79) 



where a and (3 are imprecisely known parameters of order 
unity. Minimizing the total error for fixed ft and 7 im- 
plies that the filters should be chosen to have a linewidth 
scaling as 



2r ~ o 2/3 7 1/3 



(4.80) 



where this expression would be exact for a = j3. This 
optimal scaling is interesting in that it differs from the 
geometric mean ft 1 /^ 1 / 2 which is what one m ight have 
guessed. Substituting this back into Eq. ( 4.79 ) gives 



ftotal 



2/3 



(4.81) 



where this expression would be exact for a = = 1. 
Thus with ft = 7OO7, which is readily achievable experi- 
mentally, the atomic dynamics would agree with those of 
the dressed state model with an accuracy of about 99%. 

Accurately testing these results via stochastic simu- 
lations with reasonable computational resources would 
require a long time, or highly specialized algorithms and 
extensive coding. However, the scaling laws can be tested 
in the following approximate way. From Fig. 4 we see 
that the state after a sideband detection alters little until 
the next detection into a different sideband. Therefore if 
we calculate the average state of the atom after a partic- 
ular sideband detection we should get a decent estimate 
of how close the atom usually is to a dressed state. The 
way this average state can be calculated is presented in 
Appendix A. Let us consider a high-frequency sideband 
detection (through cavity a), and call this average state 
p a - We expect this state to be close to the low-energy 
dressed state |— )(— |. Therefore we can define the ap- 
proximate probability of error as 



c app 



Pa 



(4.82) 



An expression for this probability for error is derived 
in Appendix A, and a perturbation expansion in 7/T and 
T/f] yields 



sr 2 

4ft 2 



J_ 
8T' 



(4.83) 



Comparing to expression ( |4.79| ) above shows that we have 
a = 5/4 and j3 = 1/2, which are of order unity as ex- 
pected. Since the two methods for calculating the proba- 
bility of error give the sa me sca ling, we c an be confident 
in the final results of Eq. ( [4.80D and Eq. ( |4.8l| ). 



V. TWO-STATE JUMPING 

We have found that the dressed atom theory rather 
than the TM theory, well approximates the evolution of 
the atom under perfect spectral detection in the high- 
driving limit ft 7. However, this does not prove that 
the TM theory does not does not describe the atomic dy- 
namics under some other detection scheme which Teich 
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and Mahler failed to identify. To investigate this ques- 
tion we turn now to the se cond of the three features of 
the TM model listed in Sec. II B. That is, in steady state, 



the atomic state is always one of two fixed pure states. 



A. Homodyne Detection 

Since we are seeking a physical detection scheme which 
will yield two-state jumps, the theory of quantum trajec- 
tories expounded in Sec. IV A is again the appropriate 
theory. For the system to remain a pure state it is neces- 
sary for the detection scheme not to entangle the atomic 
state with any other systems (as occurs in spectral de- 
tection). Thus we do not require cascaded quantum sys- 
tems theory and instead consider measurement operators 
fl a (dt) in the Hilbert space of the atom alone. 

With this restriction it might be thought that the only 
option is then the measurement operators 



Qo{dt) = 1 — i—a x dt ~ —a^adt, 
Qi(dt) = \J jdt a. 



(5.1) 
(5.2) 



discussed in Sec. IV A, which correspond to direct detec- 
tion and give rise to quantum jumps quite unlike those of 
the TM theory. However this is not the case. Although 
the master equation ( |2.lD i s invar iant under the transfor- 
mation given by E ns. ( | 2.16| ) a nd (2.17 ), the measurement 
operators of Eqs. (4.9|) and ( 4.1C ) are not invariant. In 



the context of the two-level atom, the master equation 



■ n i 1 



(5.3) 



can be unraveled by any of a family of measurement op- 
erators parameterized by the complex number /i, 



tt (dt) = 1 - I i—a x 
ni(di) = yj^fdt (a + fi). 



7 f . * 



71/4 



dt, (5.4) 
(5.5) 



Direct detection is recovered by setting fj, = 0. 

The transformation parameterized by [i can be physi- 
cally performed simply by adding a coherent amplitude 
to the field radiated by the atom before detecting it. If 
the atom radiates into a beam as considered previously, 
this can be achieved by mixing it with a resonant local 
oscillator at a beam splitter, as shown in Fig. 5. This is 
known as homodyne detection. The transmittance of the 
beam splitter must be close to unity and the local oscil- 
lator strength chosen such that the transmitted field (in 
the absence of the atom) would have a photon flux equal 
to 7|/i| 2 - The phase of /i is of course defined relative to 
the field driving the atom. 

Since our aim is for the atom to remain in one of two 
fixed pure states, except when it jumps, we must examine 
the fixed points (i.e. eigenstates) of the operator fl (dt). 



It turns out that it has two fixed states, such that if 
Re[/x] 7^ 0, one is stable and one unstable. For Re[/i] > 0, 
the stable fixed state is 



\r'J) = ( \/0 2 -2i07 M *-^ + ^ |<?> + fi|e). (5.6) 



Here the tilde denotes an unnormalized state. The cor- 
responding eigenvalue is 



K = -7 



2 n 2 _ 2l n 1 p* - ^. (5.7) 



The unstable state and eigenvalue are found by replacing 
the square root by its negative. 



B. Adaptive Homodyne Detection 

Let us say jjL = with Re[/i + ] > 0, and assume the 
system is in the appropriate stable state When a 

jump occurs the new state of the system is proportional 
to 



n+|^)«( £ r + A* + )K+>. 



(5.8) 



The new state will obviously be different from ) and 
so will not remain fixed. This is in contrast to what we 
are seeking, namely a system which will remain fixed be- 
tween jumps. However, let us imagine that immediately 
following the detection, the value of the local oscillator 
amplitude /i is changed to some new value, This is 
an example of an adaptive measurement scheme | p8|,p9[ , 
in that the parameters defining the measurement depend 
upon the past measurement record. We want this new 
to be chosen such that the state (a + n + )\ip+) is a 
stable fixed point of the new ^(dt). The conditions for 
this to be so will be examined later. If they are satis- 
fied then the state will remain fixed until another jump 
occurs. This time the new state will be proportional to 

(a + H_)(a + = [fx-fi- + (fi- + n+)a]\ip+). 

(5.9) 

If we want jumps between just two states then we require 
this to be proportional to Clearly this will be so if 

and only if 



-fi + . 



(5.10) 



Writing /i+ = /i, we now return to the condition that 
((T+fi + )\ip+) be an the stable fixed state of f2j7(dt). From 
Eq. (5.6), and using Eq. ( 5. 10| ) , this gives the relation 



n 2 + 2mjp,* - -j = y n 2 - 2^ 7/ i* - ^- — . 



(5.11) 



This has just two solutions, 
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(5.12) 



which, remarkably, are independent of the ratio 
The stable and unstable fixed states for this choice are 



\g)- 



and the corresponding eigenvalues are 



7 in 

'8 2' 
57 iO 

'T T T' 



Note that the unstable states \tp^_ 
dressed states |±) defined in Eq. (3. 
f2 3> 7, the stable states l?/^ 
nal dressed states |=F 



= |e), (5.13) 
(5.14) 

(5.15) 

(5.16) 

are simply the 
In the limit 



become equal to orthogo- 
Like the states \h), \ l) defined in 
Eqs. (4.24), (4.25), they differ from the dressed states by 
an amount of order (7/fi) 2 . However they are different 
states, as their locus on the Bloch sphere in Fig. 6 shows. 
In any case, the system evolution in steady state, jump- 
ing between the stable states, corresponds closely to the 
dressed state evolution. It is shown in Appendix B that 
the system rapidly reaches this stationary evolution from 
an arbitrary initial condition. 

Unlike in spectral detection, the atomic state is not 
entangled with any other system and does jump cleanly 
from one pure state to another. The total error, the prob- 
ability for the atom not to be in the nearest dressed state, 
is just 



Etotal 



7 



Ail 2 + 2 7 2 



(5.17) 



which goes as (7/2^) 2 in the limit of strong driving. Note 
that as 7/fi goes to zero, this error goes to zero much 
faster (with power 2 compared to power 2/3) than the 
corresponding minimum error for spectral detection in 
Eq. ( pL79l ). 

This adaptive measurement scheme described above 
would be relatively easy to implement experimentally 
(assuming that the problem of collecting all of the atomic 
fluorescence has been solved.) It requires simply an am- 
plitude inversion of the local oscillator after each detec- 
tion. The other difference from usual homodyne detec- 
tion is that the transmitted local oscillator intensity is 
very small: it corresponds to half the photon flux of the 
atom's fluorescence if the atom were saturated. In ei- 
ther stable fixed state, the actual photon flux entering 
the detector in this scheme is 

(tpt I (M± + + "M*) = \, (5-18) 

which is again independent of il/7- This rate is also of 
course the rate for the system to jump to the other sta- 
ble fixed state. The rate 7/4 is the same as the rate 



of state-changing transitions in the dressed atom model. 
However, since there are no detections which leave the 
atomic state unchanged, the total rate of photedetections 
is half that of the dressed atom or TM model. 



VI. ORTHOGONAL JUMPING 

The final detection scheme we will examine in this pa- 
per is quite similar to the two-state detection scheme. 
Rather than reproducing the two-state feature of the 
TM model, it aim s to reproduce instead the third fea- 
ture listed in Sec. |ll B| . That is, it will ensure that at 
all times, the state after a jump is orthogonal to the one 
before. This can be achieved using homodyne detection 
as in the two-state model. The condition is 



<V>|(a + M M = o, 



/i(t) = -(V(t)klV>(*)>- 



(6.1) 



(6.2) 



Like two-state jumping, this orthogonal jumping 
clearly requires an adaptive measurement scheme in that 
the amplitude (in this case intensity and phase) of the 
local oscillator will depend on the previous measurement 
history, which will determine the current state \ip(t)). 
But since in this case the amplitude \i is found directly 
from \ip(t)), it will also depend on the initial state of the 
system at some time in the past |^(0)). That is, it re- 
quires the experimenter to know the initial state of the 
system. This is an unusual condition that will be dis- 
cussed more later. 



A. Dynamics 

The "jump" measurement operator corresponding to 
the choice in Eq. (|6.2|) is 



fi x (eft) = ^dt [a ~ (a)(t)] 



(6.3) 



From Eq. (5.4), the no-jump measurement operator is 

7lM| 2 



n (di) = 1 



■ n . 7 + , / \* 
1 — u x H — er' a + 7(07 a 



dt. 
(6.4) 



For f2 > 7, this nonlinear operator has the following sta- 
ble fixed states: 



V2n\e ± ) = n\ g ) + (i 7 ±^ 2 -7 2 ) k) 



(6.5) 



Once again in the limit that Q/7 becomes large, these ap- 
proximate the dressed states, differing from them by an 
error of order /2Vt) 2 . However, as shown in Fig. 6, they 
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are different both from the states \h), \l) appearing in our 
analysis of spectral detection and the two stable states 
in the two-state jumps. Linearizing the nonlinear 
operator fio(dt) about either of the fixed states gives the 
eigenvalues 

A = -J±^/(£) 2 + (7 2 -tn (6.6) 

For CI > -y/17/167, the eigenvalues are complex, so that 
the two fixed points on the Bloch sphere are foci of the 
no-jump dynamics. For SI > 7, the real and imaginary 
parts of the eigenvalues approach —7/4 and Q respec- 
tively. 

Because of the nonlinearity of the evolution described 
by the above Cli(dt) and £lo(dt), an analytical approach 
is more difficult in this case. Indeed, for f2 ~ 7 we find 
numerically that the evolution is very complicated. How- 
ever, for fl ^> 7, the two stable fixed states approach or- 
thogonality. This means that if the system is in a stable 
fixed state, then a jump (which by construction takes the 
system to an orthogonal state) will take it close to the 
other stable fixed state. Thus we expect the two-state 
jumping of the previous section to be approximately re- 
produced in the strong-driving limit. This expectation is 
confirmed by the numerical simulations shown in Fig. 7. 
As expected from Eq. ( |6.6D , the state after a jump spirals 
towards the closest fixed state. 

Once again, attempting to replicate a feature of the 
TM model has in fact replicated approximately the be- 
haviour of the simple dressed atom model. The devia- 
tion of the orthogonal jump behaviour from that sim- 
ple model can be roughly estimated as follows. First, 
as noted above, the stable fixed points differ from the 
dressed states by an error of order (7/2SI) 2 . Second, 
when a jump from \0+) occurs the new state is such that 
the dressed state |— ) lies midway between it and \0~). 
The error immediately after such a jump is thus also 
(j/2fl) 2 . Hence we can estimate that the overall error 
(the probability of the atom not being in the expected 
dressed state at steady state) is 

From numerical simulations shown in Fig. 8 we find 




so that the error is roughly three times greater than in 
the in the two-jump case. 

B. Nonlinearity and Consistency 

The "orthogonal jump" evolution presented here has 
been considered before. Breslin et al. Q used it to ex- 
amine questions of information production and quantum 



chaos. Earlier, it was suggested by Diosi H] as a unique 
way of unraveling a master equation. 

The orthogonal jump evolution differs from the other 
unravelings considered here in that it is nonlinear in the 
sense that the measurement operators Qi(eft), Qo(dt) de- 
pend upon the state of the system. This is only possible 
if we assume that the experimenter who is monitoring 
the system knows what its initial state is. Another ex- 
perimenter, arriving in the middle of the monitoring and 
having not communicated with the first experimenter, 
would at assign a different (mixed) state to the system. 
That second experimenter would then disagree with the 
first experimenter's control of the measurement appara- 
tus (the local oscillator amplitude) because the two ex- 
perimenters would assign different values of (c) to the 
system. 

This disagreement between two observers on how the 
system will (or rather, should) evolve applies also the TM 
model. This was pointed out in Ref. p8[ , where it was 
argued that this was a fatal flaw in the consistency of the 
TM model. In the present context we can now see that 
the argument in Ref. [ fl8| is not wholly convincing. It is 
possible for different observers to disagree on the state of 
the system, but as long as one observer "holds the reigns" 
of the equipment, the future evolution is unambiguous. 

Thus, while the TM model fails on other grounds (in 
that its cannot be physically realized), its internal con- 
sistency is technically no worse than that of the orthog- 
onal jump model. However, Teich and Mahler originally 
claimed that their model is "the stochastic process which 
governs the time evolution of an individual quantum sys- 
tem" . That is, it is apparently supposed to represent an 
observer-independent reality. In this spirit, the disagree- 
ment of two observers on the dynamics of the system is 
still a serious problem. 

VII. CONCLUSION 

In this work we have compared the evolution of a 
strongly driven, damped two level atom under a variety of 
stochastic evolution models. The model which served as 
the starting point for our investigations was the one due 
to Teich and Mahler. In the steady state, this predicts 
that the atom jumps between two orthogonal states, 

s/2\4> l ) = \ 9 )-i\e), (7.1) 
V2\<f> 2 ) = \g)+i\e), (7.2) 

with rates approximately equal to 7/4. Here 7 is the 
spontaneous emission rate for the atom. The states 
\4>i)i |02) are the states which diagonalize the stationary 
state matrix for the system in the limit of strong driving. 

The TM model is supposed to be an objective descrip- 
tion of the behaviour of individual quantum systems. 
However, this objectivity runs counter to one of the fun- 
damental features of quantum mechanics, entanglement. 
A fluorescent atom becomes entangled with the state of 
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the field into which it emits, and different ways of moni- 
toring that field will give different information about the 
atom and hence collapse the atom into different states. 
These different processes are called unravelings of the 
master equation of the atom. The question we then posed 
is, can any unraveling reproduce the theory of Teich and 
Mahler? 

The short answer is no. The long answer is far more 
interesting. Individual features of the TM model can 
be mimicked. First, Teich and Mahler claimed that the 
jumps of their atom corresponded to the emission of pho- 
tons with different frequencies. We modeled this process 
using the dressed atom model (Sec. Ill) and then a more 
exact and far more complicated method explicitly includ- 
ing the filters used for distinguishing the different fre- 
quencies (Sec. IV). The rate and ordering of jumps were 
roughly as predicted by the TM model. Then (Sec. V), 
we derived the measurement scheme which ensures that 
the atom, in steady state, jumps between precisely two 
states (as in the TM model). Last (Sec. VI), we derived 
the scheme which ensures that when the atom jumps to 
a new state, it is orthogonal to the old one (as in the TM 
model). 

What is interesting is that in all of these cases the atom 
spends most of its time close to one of the following two 
states 



V2\±) = \g)±\e). 



(7.3) 



These we have called, in a slight abuse of terminology, 
the dressed states of the atom. That is because in the 
simple dressed-state theory, these are precisely the states 
the atom jumps between for a coherent driving field. The 
probability for error, that is the probability for the sys- 
tem to be in a state other than the expected one of these 
dressed states, depends on the monitoring scheme. For 
the three schemes considered here, the error probabilities 
are respectively 



3 / 7 \2/3 

^spectral ^ ^ 



l n 



^two— state 



^orthogonal ~ ^ ^ ^ 



4 

3 /7^ 2 



(7.4) 
(7.5) 
(7.6) 



where f2 3> 7 is the Rabi frequency of the driving field. 

These results are in flat contradiction to those of the 
TM model. The dressed states are as far away as possi- 
ble from the diagonal states \4>i), \4>2)- But our purpose 
here is not to bludgeon the TM model by repeated in- 
stances of nonphysicality, but to marvel at the fact that 
whenever one tries to mimic it, one ends up instead mim- 
icking the evolution of the dressed state model. It seems 
as if jumping between the dressed states is the evolution 
which the atom "wants to do" . Of course one can force 
it to behave otherwise. Direct photodetection will result 
in quite different evolution. But attempting to make the 



evolution simple seems to lead inevitably to the dressed 
state jumping. 

The lesson here is that diagonal states, that is, the 
states which diagonalize the state matrix, have no rela- 
tion to the states the system prefers. In the case of a res- 
onantly driven two-level atom the preferred states seem 
to be the dressed states. Lest the reader be annoyed at 
our repeated attribution of state preference to a quan- 
tum system, we point out that this terminology is widely 
used in the study of decoherence and the classical limit 
p2| . A more technical term is "environmentally- induced 
superselection" , which emphasizes that the preference of 
certain quantum states is a property of the environment 
of the system, as well as the system itself. 

A two-level atom is such a small quantum system that 
it might be thought that there is no point considering a 
classical limit. This is a fair comment if one supposes 
that a classical limit must be a deterministic limit. But 
if one is prepared to accept a stochastic classical model, 
than it seems that something like the dressed state model 
is the appropriate limit for a strongly driven atom. 

Of the measurement schemes analyzed here, the two- 
state scheme of Sec. V gives the best approximation to 
the evolution arising from the dressed state model. A rig- 
orous way of defining the preferred unraveling of a mas- 
ter equation has been formulated recently by one of us 
P3| , using the concept of robustness. Preliminary work 
suggests that the two-state jumping scheme is in fact 
the most robust unravelling of the resonance fluorescence 
master equation. This will be pursued in a future work. 
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APPENDIX A: AVERAGE STATE 
CONDITIONED ON A FILTERED PHOTON 
DETECTION 

In this appendix we show how to calculate the aver- 
age state of the atom given that a detection through a 
particular filter has just occurred. This calculation re- 
quires consideration of only a single filter cavity, which 
is described by the master equation 



W 



-i[(Q/2)a x ,W]+jV[(r}W 
-i [u a <Ja,W] +2TV[a]W 
- VF7 (oTTV + aWa) - a)aW 



Waa 



(Al) 



which is just a rearrangement of Eq. (4.17). The first 
line is the atomic dynamics, the second line is the cav- 
ity's free dynamics and the third line is the coupling of 
the atom to the cavity. 

The average state of the atom immediately after the 
detection of a photon passed by the cavity is 
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Pa 



Tr cav [aWa+] 
Tr [aVFa+] ' 



(A2) 



where here W is the stationary solution of Eq. flAip . To 
determine this state, we first trace out the field entirely 
to obtain the familiar master equation for the TLA of, 

^ = -i m/2)a x , po] + jV[a] Po = C Po . (A3) 

Writing the atomic state matrix in the Bloch representa- 
tion, 



Pa = 2 (Po 1 + x o&x + VaVy + z <J z ) , 
the relaxation superoperator takes the form 



(A4) 



£ = 



I 





V-7 







\ 





o -h -n 



(A5) 



n -7 J 



acting on the vector (p , Xq, yo, z ) T . Assuming that the 
state is initially normalized (po = 1), the stationary so- 
lution is 



(po,x ,yo,z ) 



2^7 

1 

' ' 7 2 + 2r! 2 '7 2 



■7 



2f7 5 



(A6) 



Next, we calculate the steady-state value of pi 
Tr cav [aW]. From Eq. <§K}§ we get 



dpi 
dt 



Cpi - iuJaPi - I>i - VT7 



ap . 



(A7) 



Converting the operator apo into the Bloch vector rep- 
resentation, we find that the steady-state representation 
of pi is 



Pi 



/T7 



(c - iiv a - r) 1 



x (x - iyo,Po + z a , -ip - iz ,-x + iy ) (A8) 
We can now repeat this process for p2 = Tr ca _ v [aWa']. 



Tracing over Eq. 



Al 



gives 



dp2 r or , 
— = C P2 2T P2 



(picrt + ap Q . ( A9 ) 
Thus the value of p 2 in terms of p\ is given by 



solution in the limit 7/r, T/f2 <C 1. We are particularly 
interested then in the case where the cavity is tuned to 
a sideband, with u> a = f2. Writing pi in the Bloch repre- 
sentation, the probability of error (which is the quantity 
we ultimately require), can be written as 



(+1/02 1+) _ X2+P2 



c app 



Tr[p 2 



2p 2 



(All) 



Then, using the symbolic series expansion capability of 
MATHEMATICA, we find 



sr 



5T 2 
4^ 2 



7 



32r 2 



7r 7 
8S12 



(A12) 



The two leading terms in this expansion are quoted in 
the main text. 

Not relevant to the problem at hand, but nevertheless 
interesting, is the question of what would happen in a 
different limit, namely T <ti 7 <C fi. Here the cavity is 
so narrow that any photon it passes is definitely from 
one peak of the Mollow triplet only (assuming the cavity 
is tuned appropriately). It might be thought that this 
would condition the state of the atom very well In fact 
the perturbative result obtained for uj a = Q + A, with 
A <C Omega, is 



P2 oc 



7 



A 2 



7)' 



4r 

1 <?x 

7 



(A13) 



The Lorentzian lineshape here is as expected from the 
analysis of Mollow ||, but the linear dependence of X2 
upon r/7 shows that the narrower the cavity the worse 
the conditioning of the state. This can be understood as 
follows. The detection of a photon having passed through 
the cavity actually gives information about the atom at 
the time the photon entered the cavity. The transit time 
of the photon through the cavity has an exponential dis- 
tribution 2Te~ 2rt . If the atom did jump into the negative 
a x eigenstate at time 0, then by time t the expectation 
value of a x would be equal to — e -7 '/ 2 . The average value 
of u x at the time of detection would thus be equal to 



4r 



2T e - 2Vt -^l 2 dt : 
7 



(A14) 



for r <C 7. This agrees with the result in Eq. ( A13 ) once 
it is normalized. 



P2 



' [x r 1 -y[,p\ + z[,-p[- z\,-x r 1 +y[) T . 



C~2T 



(A10) 



where the r and i superscripts indicate real and imagi- 
nary parts. Substituting in the above expression for p\ 
and normalizing thus gives p a as required. 

Due to the algebraic complexity of the exact matrix in- 
verses, it is more instructive to examine the perturbative 



APPENDIX B: TRANSIENT BEHAVIOUR OF 
THE TWO-JUMP MODEL 



From Sec. VB it is evident that an atom prepared in 
one of the stable states |^) will always jump between 
the two stable states. However it remains to be shown 
that the atom will relax towards this behaviour for any 
initial state. This can be done as follows. Consider 
the local oscillator set to p + and the atom in a state 
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\ip(0)) = p\4>t) + Q'lV'u")- If the system then evolves for 
a time £ before making a jump, its unnormalized state 
\ip(t)) after the jump is 



pexp(A+t) 2qexp(A+t) 



2 20 + ^7 

(2n-H)gexp(A+t) 

2(20 + i 7 ) IVu; ' 



(Bl) 



Here (ip(t)\ip(t)) equals the probability per unit time for 
the jump to occur. 

Two-state jumping will be stable if the average jump, 
at time t, causes KV'WIVv")! 2 to decrease. Conveniently 
there is no need to examine the case for — fj, because of 
the symmetry of the situation. The quantity we wish to 
calculate is the average of KV'WIV'u )| 2 ' weighted by the 
probability of the jump occurring at time t. The result 



E[|W*)|V£>| S 



\(mWu)\ 2 w)w))dt (B2) 



\(4,m-)\ 2 dt. 



(B3) 



For f2 ^> 7, the stable and unstable states are nearly 
orthogonal so we have 



E[|W*)|VV 



\qeMKt)\ 2 dt 



f- (B4) 



On average, then, after each jump the probability for the 
atom to be in the unstable state is reduced to a fifth of 
its size prior to the jump. Thus we can conclude that the 
two-state evolution is stable. 
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FIG. 1. Diagram of the experimental configuration for ob- 
taining the fluorescence of an atom as a beam, and passing it 
through a filter cavity. 



FIG. 2. Probability for the atom to be in the dressed state 
|— ) for a typical observation record from spectral detection 
using a single cavity tuned to the resonant frequency of the 
atom. The symbol a denotes the state following the detection 
of a photon at photodetector a and I at photodetector 1 (see 
Fig. 1). The parameters used are fi/7 = 50, r/7 = 8. 




Time (y ') 

FIG. 4. Probability for the atom to be in the dressed state 
|— ) for spectral detection using two cavities tuned to the side- 
bands of the Mollow triplet, ujq ± fi. The symbol a denotes 
the state following the detection of a photon at photodetector 
a, b at photodetector b, and I at photodetector 1 (see Fig. 3). 
The parameters used are fi/7 = 50, r/7 = 8. 
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FIG. 5. Diagram of the experimental configuration for a 
homodyne measurement of the fluorescence of an atom. The 
amplitude of the local oscillator is assumed to be variable as 
a function of time, determined by an electro-optic modulator 
(EOM). The modulator can be controlled by the experimenter 
using the results of the measurement. 
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FIG. 3. Diagram of the experimental configuration for 
passing the fluorescence of an atom through two filter cav- 
ities. 



FIG. 6. Locus of various states of the Bloch sphere: \h), \l) 
from Sec. IV for fi/7 £ [1/2, 00) (solid line and circle); {ip^ 1 ) 
from Sec. V for fi/7 £ [0, 00) (dashed line and square); \9±) 
from Sec. VI for fi/7 6 [1, 00) (solid line and cross). The cir- 
cle, square and cross show the states for fi = 27. Also shown 
are the dressed states |±) (asterisks) and the high-driving 
limit of the TM states \<f>i), |</>2) (diamonds). An equal area 
projection of the sphere in terms of the Euler angles <j>, 9 is 
used. 
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FIG. 7. Stationary probability distribution of the atomic 
state on the Bloch sphere under the orthogonal jump detec- 
tion scheme. The distribution is approximated by 1000 points, 
which are separated in time by 0.015 7 _1 . The ratio of driv- 
ing to damping is fi/7 = 10. The same equal area projection 
of the sphere as in Fig. 6 is used, but note that the <f> axis is 
cut at a different point. 




FIG. 8. Mean error (ie. probability for the system not to 
be in the predicted dressed state) versus fi/7 for orthogonal 
jump evolution. The fitted curve is e ap p = 2.99 x ( r y/2Q.) 2 
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